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Abstract 

In  this  paper,  an  analytical  model  of  a  micro  solid  oxide  fuel  cell  (SOFC)  system  fed  by  butane  is  introduced  and  analyzed  in  order  to  optimize 
its  exergetic  efficiency.  The  micro  SOFC  system  is  equipped  with  a  partial  oxidation  (POX)  reformer,  a  vaporizer,  two  pre-heaters,  and  a  post¬ 
combustor.  A  one-dimensional  (ID)  polarization  model  of  the  SOFC  is  used  to  examine  the  effects  of  concentration  overpotentials,  activation 
overpotentials,  and  ohmic  resistances  on  cell  performance.  This  ID  polarization  model  is  extended  in  this  study  to  a  two-dimensional  (2D)  fuel 
cell  model  considering  convective  mass  and  heat  transport  along  the  fuel  cell  channel  and  from  the  fuel  cell  to  the  environment.  The  influence  of 
significant  operational  parameters  on  the  exergetic  efficiency  of  the  micro  SOFC  system  is  discussed.  The  present  study  shows  the  importance  of 
an  exergy  analysis  of  the  fuel  cell  as  part  of  an  entire  thermodynamic  system  (transportable  micropowerplant)  generating  electric  power. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Conventional  power  generation  systems  transform  chemi¬ 
cal  energy  of  a  fuel  into  useful  electrical  or  mechanical  power 
with  an  intermediate  step  of  heat  production.  Fuel  cells  are  a 
remarkable  alternative  for  power  generation  transforming  chem¬ 
ical  energy  directly  into  electric  power  and  achieving  poten¬ 
tially  higher  efficiencies  than  conventional  power  generation 
systems.  A  promising  application  of  fuel  cells  are  micro  fuel 
cell  systems  generating  electric  power  at  the  order  of  some 
watts  to  power  mobile  devices,  including  cell  phones,  cam¬ 
eras,  and  notebooks.  Micro  fuel  cell  systems  have  the  poten¬ 
tial  to  replace  rechargeable  batteries  used  in  mobile  devices 
today. 

A  traditional  method  of  analyzing  power  generation  systems 
is  an  energetic  analysis  applying  the  first  law  of  thermodynam¬ 
ics.  However,  it  is  clear  that  instead,  an  exergetic  analysis  with 
exergy  as  the  measure  of  the  quality  (useful  part,  transformable 
to  work)  of  energy  can  be  used  to  specify  design  optima  which 
are  different  from  those  resulting  from  the  energy  conservation 
law  [1].  SOFC  systems  have  been  analyzed  in  the  sense  of  exer- 
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getic  efficiency  over  the  last  few  years.  SOFC  systems  fed  by 
ethanol  [2-5],  methane  [4-7],  and  hydrogen  [7]  involving  differ¬ 
ent  forms  of  fuel  reforming,  i.e.  steam  reforming,  reforming  with 
carbon  dioxide,  and  POX  reforming,  have  been  investigated.  A 
point  that  needs  improvement  in  these  otherwise  interesting  stud¬ 
ies  [2-7]  is  the  simplistic  fuel  cell  model.  Improved  ID  fuel  cell 
models  accounting  for  polarization  resistances  more  accurately 
have  been  developed  [8,9]. 

The  objective  of  the  present  study  is,  first,  the  extension  of 
existing  ID  polarization  models  [8,9]  to  a  2D  SOFC  model.  This 
extended  fuel  cell  model  is  then  integrated  into  an  entire  SOFC 
system,  which  is  analyzed  using  the  approach  of  exergetic  effi¬ 
ciency  as  also  presented  in  [2-7] .  The  operational  parameters  of 
the  system  are  chosen  to  match  experimental  data  from  previous 
studies.  The  objective  of  this  study  is  the  analytical  investiga¬ 
tion  of  the  exergetic  efficiency  of  the  system  and  it  is  beyond  the 
purpose  of  the  present  paper  to  compare  the  model  with  a  real 
operating  SOFC  system.  Specific  requirements  of  an  application 
such  as  a  micro  SOFC  system  generating  electric  power  of  some 
watts  are  considered  in  this  study.  The  developed  system  is  fed 
by  butane  and  involves  POX  reforming.  Butane  is  chosen  as  fuel 
because  of  its  relatively  easy  and  efficient  storage  in  liquid  phase 
like  in  small  cartridges.  Since  the  system  has  a  modular  design, 
different  kinds  of  reformer  as  well  as  different  fuels  could  be 
easily  integrated.  This  study  puts  forth  possibilities  to  increase 
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Nomenclature 

a  flow  availability  (W) 

a  molar  availability  (J  mol  - 1 ) 

^chem  molar  chemical  availability  (J  mol-1) 

A  exergy  flow  rate  (W) 

b  width  (m) 

c  molar  concentration  (mol  m-3) 

D  diffusion  coefficient  (m2  s-1) 

E  cell  voltage  (V) 

E°  Nernst  potential  (V) 

Tsa,  Ec  activation  energy  (J  mol-1) 

Ee\  activation  energy  (J  mol  - 1 ) 

F  Faraday  constant  (96485  C  mol- 1 ) 

h  absolute  enthalpy  including  enthalpy  of  formation 

(W) 

h  molar  absolute  enthalpy  including  enthalpy  of 

formation  (J  mol-1) 
i  current  density  (Am-2) 

i{)  exchange  current  density  (A  m-2) 

/q  a,  /q  c  pre-factor  of  exchange  current  density  (A  K  m-2) 
j  index 

J  molar  flow  (mol  m-2  s- 1 ) 

/  thickness  (m) 

L  cell  length  (m) 

M  molar  mass  (g  mol- 1 ) 

n  molar  flow  rate  (mol  s  - 1 ) 

nQ  number  of  electrons 

p  pressure  (Pa) 

P  power  (W) 

Pe  Peclet  number 

Pr  Prandtl  number 

Q  heat-transfer  rate  (W) 

r  average  pore  radius  (m) 

R  universal  gas  constant  (8.3145  J  mol-1  K-1) 

Rq i  area- specific  ohmic  resistance  (Q  m2) 

s  entropy  (WK-1) 

s  molar  absolute  entropy  (J  mol- 1  K- 1 ) 

S  surface  area  (m2) 

T  absolute  temperature  (K) 

u  velocity  (ms-1) 

U  utilization  factor 

v  along-the-cell  direction  (m) 

X  mole  fraction 

y  across-the-cell  direction  (m) 

z  number  of  electrons 

Greek  letters 

a  heat- transfer  coefficient  (W  m-2  K- 1 ) 

8  ratio  of  diffusion  coefficients 

£  extension  of  reforming 

y  porosity 

t]  overpotential  (V) 

cpco  molar  ratio  between  CO  and  CO2 
(pn2  selectivity 


k  thermal  diffusivity  (m2  s-1) 

X  molar  ratio  of  inlet  oxygen  and  fuel 

/i  exergetic  efficiency 

o 

a  collision  diameter  (A) 

cro ,  (Tei  ionic  conductivity  (S  m- 1 ) 

§  tortuosity 

Q  collision  integral 

Subscripts 

a  anode 

act  activation 

c  cathode 

cone  concentration 

dest  destruction 

eff  effective 

el  electrolyte 

h  channel 

heat  heat  transfer 

in  inlet 

irr  irreversibility 

K  Knudsen 

loss  loss 

ohm  ohmic 

out  outlet 

pc  post-combustor 

pox  partial  oxidation  reformer 

pre  pre-heater 

prod  production 

rel  released 

s  surface 

sofc  solid  oxide  fuel  cell 

vap  vaporizer 

0  at  the  standard  reference  state 

(p=  1  atm,  To  =  298  K) 

Superscripts 
a  anode 

c  cathode 

sat  saturation 

mole  specific 


the  exergetic  efficiency  of  a  micro  SOFC  system  by  optimizing 
the  operating  conditions. 

A  major  issue  is  to  show  the  importance  of  analyzing  the 
entire  fuel  cell  powerplant  system  in  the  sense  of  exergetic 
efficiency  instead  of  dealing  with  the  fuel  cell  as  an  isolated 
component.  Considering  the  micro  SOFC  system  as  a  small  and 
mobile  power  generating  device,  aspects  such  as  storage,  vapor¬ 
ization,  and  reforming  of  the  fuel  as  well  as  heat  transfer  between 
components  within  a  small  space  are  essential.  Analyzing  only 
the  fuel  cell  unit  would  mean  neglecting  these  factors.  The  deter¬ 
mining  property  of  a  mobile  energy  conversion  device  is  the 
exergetic  comparison  between  stored  fuel  and  generated  electric 
power.  High  exergetic  efficiency  means  less  fuel  consumption, 
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smaller  storage  tanks,  and  less  weight  of  the  system,  leading  to 
a  marketable  design. 

2.  Mathematical  formulation 


This  section  presents  the  governing  equations  used  to  model 
the  micro  SOFC  system.  In  Section  2.1,  a  2D  fuel  cell  model  is 
developed.  Section  2.2  introduces  a  model  of  the  entire  micro 
SOFC  system. 

A  major  approach  of  this  study  is  to  investigate  the  exergetic 
efficiency  of  the  fuel  cell  system.  An  exergy  balance  for  a  con¬ 
trol  volume  at  steady  state  is  formulated  to  calculate  the  exergy 
destruction  Adest  of  a  component  [10]: 

Ades,  =  fl  -  y)  Gin  -  U  -  yj  Gout  -  P 

+  _  y^outjfloutj,  (1) 

j  j 


using  the  molar  flow  availability  a  and  the  molar  flow  rate  n  for 
all  species  j  in  the  inlet  and  outlet  gas  mixture.  Note  that  Ts  is 
the  temperature  of  the  surface  at  the  boundary  of  the  compo¬ 
nent  where  the  heat  fluxes  Q  take  place  and  Tfl  is  the  ambient 
temperature.  P  is  the  power  generated  in  the  control  volume. 
The  subscripts  “in”  and  “out”  refer  to  the  inlets  and  outlets  of 
the  system  component,  respectively.  The  definition  of  the  molar 
flow  availability  reads  [10]: 


aj  =  0 hj(T )  -  hj(T0))  ~  T0(sj(T)  -  sj(T0)) 
+  RTq  In  (Xj)  +  af  -m. 


Note  that  Xj  is  the  mole  fraction  of  species  j.  The  molar  chem¬ 
ical  availability  fichem  is  defined  in  [11].  The  absolute  molar 
enthalpies  h  and  molar  entropies  s  include  the  enthalpy  and 
entropy  of  formation,  respectively,  and  are  determined  from  the 
polynomial  equations  [7]: 


h 


RT 


P 4  3  ,  /%r4  ,  Pf\ 
4  5  T ) 


and 


s  =  R  \  Pi  ln(T)  +  p2T  +  y  T2  +  y  T3  +  y  T4  +  Pi 


(4) 


where  the  coefficients  can  be  found  in  JANAF  tables  [12]. 


2.7.  Fuel  cell 


In  this  study,  ID  across-the-cell  models  developed  by  Zhu 
and  Kee  [8]  and  Chan  et  al.  [7,9]  are  used  and  extended  to  a 
2D  along-the-channel  model.  These  existing  ID  models  [7-9] 
represent  a  unit  cell  of  a  membrane  electrode  assembly  (MEA) 
and  describe  the  diffusive  mass  transport  and  electrochemical 
processes.  Fluid  flow  and  species  transport  along  the  channels 
as  well  as  heat  transfer  are  neglected  in  these  models.  In  order 
to  get  a  better  idea  of  a  realistic  fuel  cell,  a  ID  polarization 
model  is  incorporated  as  a  submodel  into  a  2D  fuel  cell  model 
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Fig.  1.  Schematic  of  the  fuel  cell. 


considering  convective  mass  transport  along  channels  and  heat 
transport  in  both  dimensions.  Local  current  densities  can  change 
substantially  along  a  fuel  cell  channel  as  a  result  of  fuel  and 
oxygen  depletion.  ID  models  disregard  this  aspect. 

2.7.7.  Across-the-cell  direction 

The  fuel  cell  is  modeled  as  a  system  of  three  components:  a 
cathode  channel,  an  anode  channel,  and  an  MEA.  The  MEA  con¬ 
sists  of  a  cathode,  an  anode,  and  an  electrolyte  in  between.  Fuel 
and  air  are  transported  along  the  channels  (x-direction  in  Fig.  1). 
Due  to  the  high  Peclet  numbers  Pep  =  uL/D  and  PeK  =  uLhc,  mass 
and  heat  diffusion  along  the  channel  direction  can  be  neglected. 
In  this  study,  the  Peclet  number  Pep  is  at  the  order  of  5  x  103  for 
the  anode  and  cathode  channels,  and  the  thermal  Peclet  number 
PeK  is  at  the  order  of  103.  In  the  y-direction  across  the  MEA, 
reactants  and  products  of  the  fuel  cell  reactions  are  transported 
by  diffusion  (see  Fig.  1).  The  gas  mixtures  in  the  cathode  and 
anode  channels  (II  and  III)  are  assumed  to  be  well  mixed.  There¬ 
fore,  the  conditions  at  the  boundary  of  the  channels  are  the  same 
as  the  bulk  conditions.  The  temperature  only  varies  along  the  x- 
direction,  it  is  assumed  to  be  constant  in  the  y-direction.  Oxygen 
diffuses  from  the  cathode  channel  (II)  through  the  porous  cath¬ 
ode  (IB).  At  the  interface  between  the  porous  cathode  (IB)  and 
the  electrolyte  (IA),  oxygen  is  reduced  to  oxygen  ions  O2-  and 
the  ions  are  transported  through  the  electrolyte.  On  the  other 
side,  hydrogen  diffuses  from  the  anode  channel  (III)  through 
the  porous  anode  (IC)  to  the  electrolyte/anode  interface,  where 
hydrogen  reacts  with  oxygen  ions  to  form  water,  the  product  of 
the  global  fuel  cell  reaction. 

The  cell  voltage  is  mainly  affected  by  activation  barriers  for 
electrochemical  reactions  and  by  resistances  to  mass  transfer. 
Mass  transfer  includes  the  diffusion  of  gases  in  the  porous  elec¬ 
trodes  and  the  transport  of  oxygen  ions  in  the  electrolyte.  The 
resistances  are  called  activation  polarizations  of  the  anode  and 
the  cathode,  concentration  polarizations  of  the  anode  and  cath¬ 
ode,  and  ohmic  polarization  of  the  electrolyte. 

The  cell  voltage  E  can  be  calculated  using  the  energy  balance 
of  a  differential  part  of  the  fuel  cell, 

d P  +  d(2h  T-  d/z£  "f"  =  0’ 


(5) 


336 


N.  Hotz  et  al.  /  Journal  of  Power  Sources  158  (2006)  333-347 


and  the  entropy  balance, 

djprod  -  T~ldQh  -  d*h  -  d.vj|  =  0, 


the  electrolyte  r]ohm  are  modeled  based  on  [7-9].  The  detailed 
model  is  listed  in  the  appendix. 


where  d Qh  and  d P  are  the  differential  heat-transfer  rate  and 
the  differential  power,  respectively.  The  differential  dsprod  is  the 
differential  entropy  production  due  to  irreversible  processes  in 
the  fuel  cell.  The  cathode  and  anode  flows  consist  of  different 
species  j.  The  enthalpy  and  entropy  are  calculated  by: 

^  =  Y^nJj  (7) 

j 

and 

sh  =  Ynfsj,  (8) 

j 


2.1.2.  Along-the-channel  direction 

2. 1.2.1.  Membrane  electrode  assembly.  A  mass  balance  of  the 
molar  diffusion  fluxes  across  the  boundaries  of  the  ME  A  reads: 


7h2  =  ^h2o  =  2/o2, 


(16) 


using  stoichiometry  from  Eq.  (11).  The  molar  fluxes  J  are  pro¬ 
portional  to  the  current  density  /,  i.e., 


(17) 

(18) 


respectively,  with  molar  flow  rates  indicated  by  n.  The  absolute 
molar  enthalpies  h  and  molar  entropies  s  include  the  enthalpy 
and  entropy  of  formation,  respectively,  and  are  determined  from 
Eqs.  (3)  and  (4).  Combining  Eqs.  (5)  and  (6),  the  differential 
electric  power  d P  reads: 

d P  =  -(d hi  -  Tdsch)  -  (d hi  -  Tds&  -  rdsprod.  (9) 

Since  d P  is  the  electric  power  generated  by  the  fuel  cell,  it  can 
be  written  as: 


and 


i 


(19) 


The  energy  balance  of  the  MEA  takes  into  account  enthalpy 
transport  by  diffusion  and  the  generated  electric  power.  There 
is  no  heat  flux  between  the  MEA  and  the  surrounding  channels, 
since  they  are  at  the  same  temperature.  According  to  Fig.  2,  the 
energy  balance  can  be  formulated  as: 


dP  =  zFEdn±i2,  (10) 

where  z  is  the  number  of  electrons,  dnu2  Is  the  molar  flow  rate 
of  hydrogen  participating  in  the  reaction  in  the  differential  fuel 
cell  part,  and  F  is  the  Faraday  constant  [9] .  The  global  fuel  cell 
reaction  reads: 


2H2  +  02^  2H20. 


(11) 


The  differential  electric  power  is  therefore  given  as: 


d  P  =  2  FE  dnn2. 

Combination  of  Eqs.  (9)  and  (12)  results  in: 

E  =  __ L  (JK  +  +  T_  fdsy  +  dsi_\ 

2 F  \d^H2  dnn2  J  2 F  yd/iH2  d/in2  / 

T  di'prod 
2 F  d^n2 


(12) 


(13) 


The  first  two  terms  of  the  right-hand  side  represent  the  ideal  cell 
potential,  the  Nernst  potential  E°,  and  the  last  term  describes 
the  voltage  drop  due  to  polarizations  in  the  fuel  cell.  The  cell 
voltage  E  reads  therefore  as  the  difference  between  the  Nernst 
potential  E°  and  the  sum  of  all  polarizations  ritou  given  by: 

E  =  E°-  %ot.  (14) 


The  sum  of  all  polarizations  in  the  fuel  cell  rjtot  can  be  written 
as: 


Otot  —  7?act,a  T  7?act,< 


(15) 


dhJii2  +  d  hj02  =  d  hJlii0  +  dP,  (20) 

where  dh  ,  dh  jQ  ,  and  dh  /Ho0  account  for  enthalpy  transport 
by  species  diffusion,  and  dP  is  the  generated  electric  power.  Note 
that: 

d^H2  =  7H2^H2^d-^,  (21) 

d  hj02  =  Jo2ho2bdx,  (22) 

and 

d/z/H2o  =  Ju2ohu2obdx  (23) 

applies,  where  b  is  the  electrolyte  width  (measured  perpendicular 
to  v  and  y)  and  h  is  the  molar  absolute  enthalpy  in  J  mol-1 .  The 
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Fig.  2.  Schematic  of  energy  transfer  in  the  fuel  cell. 
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generated  electric  power  d P  reads: 

d  P  =  Eib  dx.  (24) 


2. 1.2.2.  Cathode  channel.  For  the  cathode  channel,  a  mass  bal¬ 
ance  and  an  energy  balance  can  be  formulated.  The  molar  flow 
rate  of  nitrogen  in  the  channel  is  constant,  since  nitrogen  is  not 
converted  in  the  cell  reaction.  The  molar  flow  rate  of  oxygen 
no2  decreases  along  the  channel  according  to: 


dnp2 

dx 


(25) 


The  differential  energy  balance  of  the  cathode  channel  can  be 
written  as  (see  Fig.  2): 


dh^  =  _dQl  _  dhjo, 

dx  dx  dx 


(26) 


2. 1.2.4.  Entire  fuel  cell.  The  heat  loss  from  the  channels  to  the 
environment  dgh  [see  Eqs.  (5)  and  (6)]  is  given  by 

dQh  =  dQch  +  dQah.  (34) 

An  exergy  balance  around  the  entire  fuel  cell,  including  the 
MEA,  the  anode  channel,  and  the  cathode  channel,  is  formu¬ 
lated  to  compute  the  exergy  destruction  Adest,  given  by: 

dAdest  _  _  L  _  To  \  d<2h  _  <±P  _  da£  _  do| 

dx  \  Tso fc(v)  J  dx  dx  dx  dx 

where  the  availability  of  the  mass  flow  in  the  cathode  and  the 
anode  channel  read: 

ai  =  (36) 

j 

and 


The  absolute  enthalpy  of  the  mass  flow  in  the  cathode  channel 
reads: 


h 


j 


(27) 


(37) 


respectively.  The  definition  of  the  molar  flow  availability  is  given 
by  Eq.  (2). 


where  j  refers  to  the  species  oxygen  and  nitrogen.  The  molar 
absolute  enthalpy  h  is  defined  by  Eq.  (3).  For  the  heat  loss  from 
the  cathode  channel  to  the  environment, 

^  =  aLhb(T,ofc(x)  ~  T0)  (28) 

applies,  where  is  the  heat-transfer  coefficient  between  the 
cathode  channel  and  the  environment. 


2. 1.2. 3.  Anode  channel.  Species  conservation  for  hydrogen  and 
water  in  the  anode  channel  leads  to: 


d^H2 

dx 


(29) 


and 


dnn2o 

dx 


=  £>7h2o, 


(30) 


respectively.  Conservation  of  energy  in  the  anode  channel  reads: 


d h\  _  _d Ql  _  d hj^  dhj^ 

dx  dx  dx  dx 


(31) 


The  absolute  enthalpy  of  the  mass  flow  in  the  anode  channel  is 
written  as: 


h 


j 


(32) 


where  j  refers  to  the  species  hydrogen,  water,  oxygen,  nitrogen, 
butane,  carbon  dioxide,  and  carbon  monoxide  in  the  anode  chan¬ 
nel.  For  the  heat  loss  from  the  anode  channel  to  the  environment, 


dQl 

-fr  =  <b(Tsofc(x)  -  To)  (33) 

applies,  where  af  is  the  heat-transfer  coefficient  between  the 
anode  channel  and  the  environment. 


2.2.  Micro  SOFC  system 

A  thermodynamic  model  of  a  micro  SOFC,  including  a  POX 
reformer,  a  post-combustor,  a  vaporizer,  and  two  pre-heaters, 
arranged  according  to  Fig.  3,  is  presented  in  this  subsection. 
The  vaporizer  is  needed  to  provide  gaseous  butane.  The  reformer 
converts  butane  by  POX  to  hydrogen  and  carbon  monoxide.  Air 
and  butane  are  preheated.  The  exhaust  gases  of  the  SOFC  can  be 
burnt  in  a  post-combustor  to  provide  heat  for  preheating  and  to 
avoid  exhaust  of  carbon  monoxide  and  unbumt  hydrocarbons. 
The  POX  reformer,  the  SOFC,  and  the  post-combustor  gener¬ 
ate  heat.  The  vaporizer  and  the  pre-heaters  require  heat,  part  of 
which  may  be  lost  to  the  environment.  The  components  of  the 
system  are  modeled  as  control  volumes,  except  the  fuel  cell, 
which  is  modeled  according  to  the  previous  subsection.  The 
heat  transport  within  the  system  is  shown  in  Fig.  3.  The  heat 
released  from  the  SOFC  is  rejected  to  the  environment,  whereas 
the  heat  released  from  the  POX  reformer  and  the  post-combustor 
is  used  to  heat  the  two  pre-heaters  and  the  vaporizer.  The  post¬ 
combustor  is  assumed  to  provide  the  heat  required  in  pre-heater  1 
and  the  vaporizer,  and  the  POX  reformer  is  assumed  to  provide 
the  heat  for  pre-heater  2.  If  the  POX  reformer  and  the  post¬ 
combustor  release  more  heat  than  required  by  the  pre-heaters 
and  the  vaporizer,  the  waste  heat  will  be  released  to  the  envi¬ 
ronment.  The  results  will  show  that  indeed  the  amount  of  heat 
released  is  higher  than  that  needed  for  the  pre-heaters  and  the 
vaporizer.  Therefore,  even  for  a  low  efficiency  of  heat  exchange 
between  the  components  of  the  micro  SOFC  system,  the  neces¬ 
sary  heat  transfer  within  the  system  can  be  achieved. 

According  to  Fig.  3,  the  heat  exchange  between  the  compo¬ 
nents  and  the  environment  can  be  formulated  as: 

Grell  —  Qpc  —  Gprel  —  Qv  ap>  (38) 

6rel2  —  Gpox  —  Gpre2>  (39) 
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Fig.  3.  Schematic  of  the  entire  micro  SOFC  system. 


and 


Grel3  —  Gprel,loss  “1“  Gpre2,loss> 


(40) 


where  the  total  heat  flux  released  from  the  system  to  the  envi¬ 
ronment  reads: 


Grel  —  Grell  T  Grel2  T  Grel3  T  Gsofcv 


(41) 


The  total  exergy  destruction  by  irreversibilities  in  the  system 
reads: 

Adest, irr  =  eldest, vap  T  Adest, prel  T  Adest, pox  T  Adest, pre2 

T  ^dest,  sofc  T  Adest, pc?  (42) 

where  the  exergy  destruction  by  heat  transfer  Adest,heat  is  calcu¬ 
lated  implicitly  by: 


Adest, heat  —  ^1  f  <34  ~\~  Cll  ~\~  Ci\2  <3\2  P  Adest, hr- 


(43) 


2.2.1.  Conservation  laws 

For  every  component  of  the  system,  a  mass  balance,  an  energy 
balance,  and  an  exergy  balance  can  be  formulated.  Conservation 
of  mass  at  steady  state  for  a  certain  system  component  requires 
that  the  inlet  flow  rates  of  hydrogen,  carbon,  oxygen,  and  nitro¬ 
gen  have  to  be  the  same  as  the  outlet  flow  rates,  i.e., 


—  y>out„ 


where  j  refers  to  the  elements  H,  C,  O,  and  N.  Other  elements  are 
neglected  in  this  study,  since  pure  butane  is  used  as  fuel  and  air  is 
assumed  to  consist  of  21%  oxygen  and  79%  nitrogen.  Humidity 
and  carbon  dioxide  in  the  air  are  neglected. 

The  pressure  is  assumed  to  be  constant  at  atmospheric  state, 
except  for  the  pressure  in  the  storage  tank  of  liquid  butane.  The 
change  of  kinetic  and  potential  energy  is  neglected.  The  energy 
balance  of  a  certain  system  component  then  reads: 


injnjnj  +  Gin  —  EWiamj  T  Gout  T  F, 


(45) 


j 


where  h  is  the  molar  absolute  enthalpy,  n  is  a  molar  flow  rate, 
G  a  heat  flux,  and  P  the  generated  power.  The  molar  absolute 
enthalpies  h  include  the  enthalpy  of  formation  and  are  calculated 
using  the  JANAF  tables  [12]  according  to  Eq.  (3). 


An  exergy  balance,  defined  by  Eq.  (1),  is  formulated  to  calcu¬ 
late  the  exergy  destruction  Adest  of  a  component  and  to  determine 
the  exergetic  efficiency  of  the  processes  in  each  system  compo¬ 
nent. 

2.2.2.  Vaporizer 

In  the  vaporizer,  the  liquid  butane  at  a  storage  pressure  p \  is 
expanded  isothermally  to  gaseous  butane  at  po.  The  vaporizer  is 
modeled  as  isothermal  at  the  ambient  temperature,  therefore  the 
temperatures: 


T\  —  T?  —  Tvap  —  To 


(46) 


are  identical.  The  molar  flow  rate  of  butane  remains  unchanged 
in  the  vaporizer,  and  the  mass  balance  of  the  vaporizer  results  in: 


^C4H10,l  =  ^C4H10,2. 

The  energy  balance  reads: 

Gvap  =  ftC4H10,l(^2  —  h\), 


(47) 


(48) 


where  h  \  is  the  molar  enthalpy  of  liquid  butane  at  pressure  p\  and 
temperature  To  [13].  The  molar  enthalpy  of  gaseous  butane  hz  is 
calculated  according  to  Eq.  (3).  The  exergy  balance  is  given  by: 


(44)  Adest, vap  =  (  1  - 


To 


T 


Gvap  +  ftQHiodO^l  -  <32), 


(49) 


vap 


where  a\  is  the  molar  flow  availability  of  liquid  butane  at 
pressure p\  and  temperature  To.  According  to  Eq.  (2),  the  molar 
flow  availability  a\  is  calculated  as: 


a\  =  On  —  h2)  -  Tq(si  -  s2 )  +  QqHio ’ 


(50) 


since  state  2  is  at  ambient  temperature  To.  The  molar  entropy 
of  liquid  butane  is  taken  from  [13]. 

2.2.3.  Pre-heaters 

The  pre-heaters  have  two  inlets  and  two  outlets  to  preheat 
two  gas  flows  without  mixing  them.  In  the  pre-heater  1,  air  and 
gaseous  butane  are  preheated  separately.  In  the  pre-heater  2,  air 
and  the  fuel  gas  mixture  are  preheated.  The  pre-heaters  allow 
heat  exchange  between  both  gas  flows.  The  temperature  at  the 
outlets  is  assumed  to  be  identical  such  that  T 3  =  T5  and  Eg  =  T9. 
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A  mean  pre-heater  temperature  is  defined  as: 

T\n  nrp  +  Tout 


Tpre  — 


in,pre 


(51) 


where  7in  pre  is  given  by  the  implicit  formulation 

^C4Hio,2^C4Hio(7in,prel)  +  n  air ,  4 h  air  (  Tin ,  pre  1 ) 

=  ^C4Hio,2^C4Hi0(^2)  +  Wair,4^air(74) 
for  pre-heater  1,  and 


(52) 


^6^6(^in,pre2)  T  ^air,7^air(2in,pre2)  —  ^/^(T^)  “1“  ^air,7^air(^7) 

(53) 


for  pre-heater  2.  The  heat  released  to  the  environment 

2 pre,loss  —  ^pre5pre(2pre  To) 


(54) 


is  assumed  to  depend  on  the  difference  between  the  mean  pre¬ 
heater  temperature  and  the  ambient  temperature,  where  apre  is 
a  heat-transfer  coefficient  between  the  surface  Spre  of  the  pre¬ 
heater  and  the  environment.  The  necessary  heat  input  for  the  first 
pre-heater  2prei  reads: 

2prel  =  2 pre  1,  loss  +  rcC4H10,3(^C4H10,3  -  ^C4H10,2) 

T  ^air,5(^air,5  —  ^air,4)  (55) 

and  the  necessary  heat  input  for  the  pre-heater  2  reads: 

2pre2  =  2pre2,loss  H"  ^  ^6,j(^8,j  —  ^6 j) 

j 

T  ^air,7(^air,9  —  ^air,7)*  (56) 

The  pre-heater  1  is  heated  by  the  post-combustor.  The  exergy 
destruction  is  given  by: 


^dest,  pre  1  —  (  1 


(2prel  2prel,loss) 


Tpre  1  / 

+  ^C4H10,2(«C4Hio,2  -^C4H10,3) 

T  ^air,5(^air,4  —  ^air,5)- 


(57) 


The  pre-heater  2  is  heated  by  the  POX  reformer.  The  exergy 
destruction  is  given  by: 


^dest,pre2  —  (  1 


(2pre2  2pre2,loss) 


7pre2  J 

T  ^  —  ^8  j)  T  ^9,air(^7,air  —  ^9,air)- 


J 


(58) 


+  5^(1  —  <^h2)H20+  ^7.pox  —  2s  j 
x  O2  +  3.76ApoxN2  +  (1—  £)C4Hio, 


2  +  (fco  5  e 

— - -tt(1  <Ph2) 

+  (pco  2 


(59) 


where  Apox  is  the  air/fuel  ratio  of  the  POX  reformer,  cpco  the 
molar  ratio  of  carbon  monoxide  and  carbon  dioxide,  <^h2  the 
selectivity  of  the  reforming  process  for  hydrogen,  and  s  the 
extension  of  reforming.  The  latter  is  the  fraction  of  butane  con¬ 
verted  in  the  reformer,  given  as: 


s  = 


^C4H10,3  -  ^C4H10,6 


^C4H10,3 


(60) 


The  ratio 


7- pox  — 


nQ25 

^C4H10,3 


(61) 


relates  the  inlet  flow  rate  of  oxygen  to  the  inlet  flow  rate  of 
butane.  The  parameter  <^co  reads: 


fco  = 


nCO,6 


(62) 


^C02,6 

at  the  outlet.  The  selectivity  of  the  reforming  process  for  hydro¬ 
gen  </?h2  is  given  by: 

(1/5)«h2  6 


<Ph2  = - 

^C4H10,3  -  ^C4Hi0,6 

A  mean  temperature  of  the  POX  reformer  is  defined  as: 

Tin  nm  +  T(5 


(63) 


Tpox  — 


m,pox 


(64) 


where  7in?pox  is  given  by  the  implicit  formulation: 

^C4Hio,3^C4Hi0(7in  ,pox  )  +  22  air,  5^  air  (Tin  ,pox) 

=  «C4Hio, 3^3(73)  +  «air, 5^5(75). 

The  energy  balance  reads: 

Gpox  +  «6^6  =  «C4H10,3^3  +  «air,5^5 


(65) 


(66) 


and  the  exergy  balance  reads: 


eldest,  pox  —  I  1 


7o 


T, 


Gpox  +  «C4H10,3«C4Hi0,3 


pox 


+  22air,5«air,5  “  ^H2,6^H2,6  ~  ^C4Hi0,6^C4Hi0,6 
-Ho2,6^02,6  -  ^N2,6^N2,6  -  ^C02,6^C02,6 
-  22(30,6*2(30,6  -  2ZH20,6^H20,6-  (67) 


2.2.4.  POX  reformer 

The  POX  reformer  has  two  separate  inlets,  one  for  preheated 
air  (state  5)  and  one  for  vaporized  butane  (state  3).  The  reforming 
reaction  is  given  as: 

C4H10  +  Apox02  T  3.76ApoxN2 

Aswro  As 

->  5^h2H2  +  - — -———CO  +  — - C02 

1  +  <^(30  1  +  <pco 


2.2.5.  Solid  oxide  fuel  cell 

For  the  SOFC,  the  model  presented  in  Section  2. 1  is  used.  The 
inlet  fuel  is  at  state  8,  and  the  inlet  air  is  at  state  9.  The  anode  and 
cathode  outlets  are  states  10  and  11,  respectively.  An  important 
performance  parameter  of  the  fuel  cell  is  the  utilization  factor 
of  hydrogen, 


22h2,8  -  ^h2,io 
22H2,8 


(68) 
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indicating  the  amount  of  hydrogen  that  is  converted  in  the  fuel 
cell.  The  parameter 


,  _  no2,9 

^sofc  — 

^H2,8/2 


(69) 


is  the  molar  ratio  between  oxygen  and  stoichiometric  oxygen 
at  the  fuel  cell  inlet.  Using  the  utilization  factor  of  hydrogen,  a 
mass  balance  for  hydrogen  in  the  anode  channel  results  in 


^h2,io  =  ^h2,8(1  -  Uh2).  (70) 

Species  conservation  for  water  in  the  anode  channel  leads  to: 


^H2O,10  =  ^H20,8  +^H2,8Uh2.  (71) 

A  mass  balance  for  oxygen  in  the  cathode  channel  requires: 

«02,9  -«H2,8£/h2 

«o2,n  = - 2 - •  (72) 

All  other  species  have  constant  molar  flow  rates  since  they  do  not 
react  in  the  fuel  cell.  The  total  heat  released  to  the  environment 
is  written  as: 


Gsofc  —  Qh  ”1“  Qh' 


(73) 


The  heat  released  from  the  cathode  and  the  anode  channel  is 
calculated  by: 


Qch  =  achb  f  (Tsofc(x)-T0)dx  (74) 

Jo 

and 

2h  =  ahb  [  (T’sofcM  -  To)  dx,  (75) 

Jo 


using  Eqs.  (28)  and  (33),  respectively.  An  essential  parameter  of 
this  study  is  the  exergetic  efficiency  fi  of  the  entire  micro  SOFC 
system,  defined  as  the  ratio  between  the  exergy  output  (i.e.,  the 
electric  power  P)  and  the  exergy  input  (i.e.,  the  flow  availability 
of  butane  at  state  1  and  air  at  the  states  4,  7,  and  12),  given  as: 


U= - - - - - - - ,  (76) 

Cl i  T  CI4  H-  (2rj  +  a\2 
where  the  electric  power  reads: 


P  =  Eb  [  i(x)  dx ,  (77) 

Jo 


according  to  Eq.  (24). 


2.2.6.  Post-combustor 

In  the  post-combustor,  the  remaining  butane,  hydrogen,  and 
carbon  monoxide  from  the  fuel  cell  exhaust  are  burnt,  such  that, 


^H2,13  =  ^CO,13  =  ^C4Hio,13  =  0. 

(78) 

As  a  result,  the  species  conservation  equations, 

^C02,13  =  4ftc4Hio,l 

(79) 

and 

^H20,13  =  5ftc4Hio,l 

(80) 

can  be  formulated  in  the  post-combustor. 

If  ^o2,4  +  ^o2,7  <  (13/2>zc4Hio,i>  the  remaining  oxygen 
from  the  fuel  cell  is  not  enough  for  total  combustion  of  butane, 
hydrogen,  and  carbon  monoxide  in  the  exhaust  gas.  In  this  case, 
more  air  is  provided  at  state  12.  The  necessary  molar  flow  rate 
at  the  inlet  12  is  given  by: 

(l  3\ 

«o2,12  =  (  y  I  «C4H10,1  -  «o2,4  -  no2j,  (81) 

such  that 


^o2;13=0-  (82) 

If  no2,4  +  ^o2,7  >  (13/2)/ic4h10,i,  the  unbumt  oxygen  from  the 
POX  reformer  and  the  fuel  cell  is  enough  for  total  combustion 
of  butane,  hydrogen,  and  carbon  monoxide  in  the  exhaust  gas. 
The  remaining  oxygen  in  the  outlet  is  in  this  case: 

«02,13  =  «02,4  +  «02,7  -  f  yj  «c4H10,1>  (83) 

where  fto2,n  =  0-  The  molar  flow  rate  of  nitrogen  in  the  outlet 
gas  reads: 


^n2,13  =  3.76(^o2,4  +  ^o2,7  +  ^o2, 12)- 
A  mean  temperature  of  the  post-combustor  is  defined  as: 

rji  _  Tin, pc  +  T13 
7pc  —  7;  ’ 


(84) 


(85) 


where  7in?pc  is  given  by  the  implicit  formulation: 

ft  10^ lo(7in,pc)  T  ri\\h  11  (7]n  pC)  T  ^7l2^12(Tin,pc) 

=  ^io^io(Tio)  +  ft\\hu(T\\)  +  ftnhniTn).  (86) 

Conservation  of  energy  leads  to: 

Gpc  +  ^13^13  =  #110*10  +  ^11*11  +  ^12^12,  (87) 

where  <2pc  is  the  heat  released  from  the  post-combustor  and  used 
to  heat  the  vaporizer  and  the  pre-heater  1 .  The  exergy  balance 
reads: 


^dest, pc  —  ji  ^  Gpc  T  ^H2,10^H2,10 

+  ^C4Hio,lodc4Hio,10  +  ^02,10^02,10  +  ^N2,10^N2,10 
+  ^C02,10^C02,10  +  ^co,io^co,io  +  ^H20,10^H20,10 
+  ^02,11^02,11  +  ^N2,11^N2,11  +  ^02,12^02,12 
+  Hn2,12^N2,12  -  ^02,13^02,13  “  ^N2, 13^N2, 13 
-  ^C02,13^C02,13  -  ^H20, 13^H20, 13 •  (88) 


3.  Numerical  solution 

3.1.  Fuel  cell 

The  fuel  cell  model  is  solved  numerically  by  first  choosing  the 
cell  voltage  E  as  an  iteration  parameter.  The  cell  voltage  must  be 
between  0  V  and  the  ideal  Nernst  potential  E°.  Since  the  result¬ 
ing  total  overpotential  rjiot  within  every  discrete  element  of  the 
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fuel  cell  depends  on  the  local  current  density  i  and  the  tem¬ 
perature  in  the  air  channel,  the  fuel  channel,  and  the  ME  A,  the 
entire  system  of  equations  needs  to  be  solved  iteratively  using 
a  Newton  algorithm.  The  integration  along  the  flow  direction 
is  accomplished  using  a  fourth-order  Runge-Kutta  method.  For 
the  calculations  presented  in  this  study,  the  results  converge  fast. 
For  a  discretization  with  10  spatial  steps,  the  relative  error  of  the 
resulting  exergetic  efficiency  is  3.5%.  For  50, 100,  and  200  steps, 
the  relative  errors  are  0.6, 0.3,  and  0.1%.  There  is  no  distinguish¬ 
able  deviation  for  discretizations  with  more  than  500  steps. 

3.2.  Micro  SOFC  system 

In  order  to  solve  for  the  entire  micro  SOFC  system,  the  fol¬ 
lowing  parameters  are  prescribed  and  assumed  to  be  constant: 
the  molar  flow  rate  of  butane  ftc4Hi0,i,  the  ambient  tempera¬ 
ture  To,  the  temperatures  of  air  and  butane  before  and  after 
pre-heating,  the  cell  voltage  E  of  the  SOFC,  the  cell  length  L 
and  width  b ,  the  extension  of  reforming  e,  and  the  oxygen/fuel 
ratios  Apox  and  Asofc.  With  the  entire  list  of  parameters  given 
in  Table  1,  the  generated  electric  power  P,  all  rejected  heat  Q , 
and  the  profiles  of  temperature,  molar  concentration,  and  current 
density  in  the  SOFC  are  numerically  calculated.  The  thermody¬ 
namic  conditions  at  all  states  are  derived. 

4.  Results  and  discussion 

In  this  study,  the  performance  of  an  SOFC  micropowerplant  is 
calculated  using  operational  parameters  based  on  experimental 
results  of  several  recent  studies,  as  shown  in  Table  1.  It  is  beyond 
the  purpose  of  this  study  to  achieve  new  experimental  results  for 
an  entire  micro  SOFC  system. 

The  parameters  of  the  MEA  are  based  on  the  parameters 
used  in  [8,9,14,15],  as  seen  in  Table  1.  For  the  anode  and  the 
cathode,  nickel-doped  lanthanum  strontium  ferrite  (FSNF)  is 
used  as  material  for  intermediate  temperatures.  The  preparation 
and  the  characterization  of  FSNF  are  investigated  in  [14].  The 
electrolyte  consists  of  a  composite  of  sodium-doped  ceria  and 
silica.  The  electrolyte  material  for  intermediate  temperatures 
and  its  preparation  are  characterized  in  [15].  The  parameters 
needed  for  matching  the  experimental  data  of  existing  SOFCs 
[16]  are  given  in  Table  1 .  These  studies  have  shown  that  satisfac¬ 
tory  performance  can  be  achieved  for  intermediate  temperatures 
between  600  and  1000  K  using  innovative  materials. 

The  pressure  is  assumed  to  be  constant  at  atmospheric  condi¬ 
tions  and  the  pressure  drop  in  the  channels  is  neglected.  Numer¬ 
ical  calculations  for  the  entire  micro  SOFC  system  for  which  a 
schematic  is  shown  in  Fig.  3  are  performed.  Fiquid  butane  at  a 
molar  flow  rate  ^c4Hi0,i  —  10-5mols-1,  ambient  temperature 
T\  =  298K,  and  pressure  p\  =  PqHio^1  =  298  K)  =  2.43  bar 
is  vaporized  isothermally  in  the  vaporizer.  Gaseous  butane  and 
air  are  heated  in  the  pre-heater  1  to  the  temperature  of  the  POX 
reformer  rpox.  The  POX  reformer  is  assumed  to  perform  isother¬ 
mally  at  Tpox  with  an  extension  of  reforming  s  =  0.9,  a  molar  ratio 
between  carbon  monoxide  and  carbon  dioxide  <^co  =  9.0,  and  a 
selectivity  for  hydrogen  cpn2  =0.9  [17].  The  SOFC  is  fed  by 
the  products  of  the  reforming  process  and  air,  where  the  molar 


Table  1 

Parameters  of  the  standard  case 


Parameter 

Value 

Reference 

Fuel 

Molar  flow  rate  of  butane  nc4Hi0,i  (f^mol  s-1) 

10 

— 

Molar  enthalpy  liquid  butane  hc4uw,i  (kJ  mol-1) 

-117.9 

[13] 

Molar  entropy  liquid  butane  Sc4Hi0,i  (1  mol-1  K-1) 

227.8 

[13] 

Molar  availability  liquid  butane  dc4Hi0,i  (kJ  mol-1) 

2805.8 

[11] 

Pre-heaters  1  and  2 

Outer  surface  area  5prc  (cm2) 

10.0 

— 

Heat-transfer  coefficient  apre  (W  m-2  K-1) 

1.0 

— 

Partial  oxidation  reformer 

Air/fuel  ratio  A,pox 

2.22 

[17] 

Extension  of  reforming  e 

0.9 

[17] 

Carbon  monoxide/carbon  dioxide  ratio  (pco 

9.0 

[17] 

Selectivity  of  hydrogen  (p^ 

0.9 

[17] 

Solid  oxide  fuel  cell 

Cell  length  L  (m) 

0.05 

— 

Cell  width  b  (m) 

0.05 

— 

Air/fuel  ratio  A,sofc 

1.0 

— 

Cell  voltage  E  (V) 

0.7 

— 

Heat-transfer  coefficient  cathode  a £  (W  m-2  K-1) 

1.0 

— 

Heat-transfer  coefficient  anode  ajj  (W  m-2  K-1) 

1.0 

— 

Membrane  electrolyte  assembly 

Average  pore  radius  r  (pan) 

0.5 

[9] 

Porosity  y 

0.3 

[9] 

Tortuosity  f 

6 

[9] 

Anode  thickness  /a  (|Jim) 

50 

[9] 

Cathode  thickness  lc  (pm) 

50 

[8] 

Electrolyte  thickness  /e i  (fxm) 

20 

[8] 

Activation  energy  of  anode  Ea  (kJ  mol-1) 

77.2 

[14] 

Pre-factor  of  anode  i(l  (1010  A  Km-2) 

15.85 

[14] 

Activation  energy  of  cathode  Ec  (kJ  mol  ) 

77.2 

[14] 

Pre-factor  of  cathode  i{){)  c  (1010  A  Km-2) 

15.85 

[14] 

Activation  energy  Ee\  (kJ  mol-1) 

72.89 

[15] 

Pre-factor  oq  (kS  m- 1 ) 

641.5 

[15] 

All  components 

Pressure  liquid  storage  tank  p\  (MPa) 

0.243 

— 

Total  pressure  p  (MPa) 

0.1 

— 

Temperature  of  environment  To  (K) 

298 

— 

Temperatures  T\,  T2,  T4,  Tj,  and  Tn  (K) 

298 

— 

Temperatures  T3,  T5,  and  To  (K) 

1000 

— 

Temperatures  7g  and  Tg  (K) 

900 

— 

Temperature  T \  q  and  T\  \  (K) 

633 

— 

Temperature  743  (K) 

800 

— 

ratio  of  hydrogen  in  the  fuel  gas  mixture  and  oxygen  in  the  air 
is  defined  by  the  air/fuel  factor  Asofc.  The  air  is  pre-heated  in  the 
pre-heater  2,  which  also  changes  the  temperature  of  the  reformed 
gas  to  assure  that  the  fuel  inlet  temperature  7s  and  the  air  inlet 
temperature  Tg  of  the  SOFC  are  the  same.  The  heat-transfer 
coefficients  between  the  environment  and  the  cathode  channel 
a £  and  the  anode  channel  a are  assumed  to  be  1  W  m-2  K-1. 
Modern  materials  and  techniques  like  vacuum  insulation  can 
achieve  small  heat-transfer  coefficients  for  high  temperatures 
around  1000  K.  However,  high-temperature  insulation  presents 
a  major  challenge  for  the  development  of  SOFC  micropower- 
plants.  The  exhaust  of  the  SOFC  is  burnt  in  a  post-combustor. 
It  is  assumed  that  heat  released  from  the  POX  reformer  and  the 
post-combustor  can  be  used  in  the  pre-heaters  and  the  vaporizer. 
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Table  2 


Energetic  results  of  the  standard  case 


Heat  transfer  to  vaporizer  <2vap 

0.214  W 

Heat  transfer  to  pre-heater  1,  <2prei 

4.58  W 

Heat  transfer  to  pre-heater  2,  <2pre 2 

1.73  W 

Heat  transfer  from  POX  reformer  <2pox 

5.31  W 

Heat  transfer  from  post-combustor  Qpc 

9.24  W 

Heat  released  from  SOFC  <2sofc 

6.80  W 

Unused  heat  from  post-combustor  Orel  1 

4.44  W 

Unused  heat  from  POX  reformer  Qvc\2 

3.58  W 

Heat  released  from  pre-heaters  <2re  13 

0.885  W 

Total  heat  released  to  environment  <2re  1 

15.7  W 

Temperature  of  vaporizer  Tvap 

298  K 

Temperature  of  pre-heater  1,  Tprei 

649  K 

Temperature  of  pre-heater  2,  Tpre2 

833  K 

Temperature  of  POX  reformer  Tpox 

1000  K 

Temperature  of  post-combustor  Tpc 

669  K 

Temperature  of  SOFC  Tsofc 

756  K 

Table  1  defines  all  relevant  operational  and  geometric  parame¬ 
ters  of  the  entire  micro  SOFC  system.  All  results  are  based  on 

these  parameters  that  represent  the  standard  case. 

Every  param- 

eter  different  from  the  standard  case  is  mentioned  in  the  text. 

Tables  2  and  3  present  energetic  and  exergetic  results  of  the 
standard  case.  In  Table  2,  the  heat-transfer  rates  within  the  sys¬ 
tem  are  quantified.  The  vaporizer  and  the  pre-heater  1  consume 
significantly  less  heat  than  the  post-combustor  releases,  and 

the  mean  temperature  of  the  post-combustor  Tpc 

is  higher  than 

the  temperatures  Tvap  and  Tprei.  The  required  heat  exchange 
between  these  components  is  therefore  possible.  The  POX 
reformer  rejects  more  heat  than  the  pre-heater  2  actually  needs, 
and  a  sufficient  temperature  difference  exists.  The  assumption 
that  the  pre-heaters  and  the  vaporizer  can  be  heated  by  the  heat 
released  from  the  POX  reformer  and  the  post-combustor  is  there¬ 
fore  feasible.  Table  3  presents  an  exergy  analysis  of  the  system. 

The  total  exergy  destruction  by  irreversibilities 

in  the  system 

Adest,irr  [see  Eq.  (41)]  is  caused  mainly  by  the  POX  reformer, 
the  SOFC,  and  the  post-combustor.  The  exergy  destruction  by 
heat  losses  Adest,heat  is  significant.  An  exergetic  efficiency  of  the 

entire  fuel  cell  system  of  19.4%  results. 

Table  3 

Exergetic  results  of  the  standard  case 

Inlet  flow  availability 

28.1  W 

Exergy  destruction  by  vaporizer  Adest,vap 

0.245  W 

Exergy  destruction  by  pre-heater  1,  Adest,prei 

0.00518  W 

Exergy  destruction  by  pre-heater  2,  Adest,pre2 

0.528  W 

Exergy  destruction  by  POX  reformer  Adest,pox 

3.29  W 

Exergy  destruction  by  SOFC  Adest,sofc 

3.75  W 

Exergy  destruction  by  post-combustor  Adest,pc 

4.51  W 

Exergy  destruction  by  irreversibility  Adest,irr 

12.3  W 

Exergy  destruction  by  heat  transfer  Adest,heat 

7.24  W 

Outlet  flow  availability  <213 

3.11  W 

Useful  electric  power  P 

5.45  W 

Exergetic  efficiency  of  the  entire  system  /z 

19.4% 

Hydrogen  utilization  factor  U h0 

100.0% 

Fig.  4.  (a)  Electric  power  P  as  a  function  of  cell  length  L  for  different  SOFC 
inlet  temperatures  T%.  Labels  indicate  T%  in  K.  All  other  operational  parameters 
are  taken  from  the  standard  case,  (b)  Exergetic  efficiency  gasa  function  of  cell 
length  L  for  different  SOFC  inlet  temperatures  T%. 

In  the  standard  case,  a  cell  length  L  =  5  cm  is  chosen.  In  order 
to  investigate  the  dependence  of  the  power  P  on  the  cell  length 
L,  the  performance  of  the  fuel  cell  system  is  calculated  for  dif¬ 
ferent  cell  lengths  L  and  various  fuel  and  air  inlet  temperatures 
T%.  In  Fig.  4(a),  the  electric  power  P  is  plotted  as  a  function 
of  the  cell  length  L  for  different  inlet  temperatures  T%.  The  star 
marks  the  result  of  the  standard  case.  The  longer  the  fuel  cell, 
the  higher  the  electric  power  P  is.  However,  the  performance 
increase  strongly  depends  on  the  fuel  cell  inlet  temperature.  For 
high  Is,  a  maximum  power  of  about  5.45  W  is  reached  for  quite 
small  cell  lengths.  A  fuel  cell  with  operational  parameters  of  the 
standard  case  could  have  a  length  of  2  cm  instead  of  5  cm  without 
any  power  loss  for  high  inlet  temperatures.  For  lower  tempera¬ 
tures  78,  the  maximum  power  of  5. 1  W  for  7s  =  850  K  and  3.0  W 
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Fig.  5.  (a)  Electric  power  P  as  a  function  of  cell  voltage  E  for  different  SOFC  inlet  temperatures  T%.  Labels  indicate  T%  in  K.  (b)  Exergetic  efficiency  gasa  function 
of  cell  voltage  E  for  different  SOFC  inlet  temperatures  T%.  (c)  Cell  voltage  E  as  a  function  of  the  current  density  i  for  T%  =  800  K.  (d)  Current  density  profiles  i  along 
the  fuel  cell  for  T%  =  800  K  and  different  cell  voltages  E. 


for  T%  =  800  K  is  reached  for  cell  lengths  above  7  cm.  Fig.  4(b) 
shows  the  corresponding  exergetic  efficiency  fi.  The  maximum 
efficiency  amounts  to  19.4%.  Even  if  7g  =  850K,  an  exergetic 
efficiency  of  18%  is  possible  for  longer  cells.  If  T%  =  800  K,  a 
maximum  efficiency  of  about  11%  results. 

Fig.  5  shows  the  effect  of  the  cell  voltage  E  on  the  perfor¬ 
mance  of  the  fuel  cell  system.  In  the  standard  case,  £  =  0.7  V 
applies.  The  electric  power  P  changes  as  a  function  of  the  cell 
voltage  £,  according  to  Fig.  5(a).  The  labels  indicate  the  SOFC 
inlet  temperature  T%.  The  star  marks  the  result  of  the  standard 
case.  The  power  P  increases  linearly  with  the  cell  voltage,  inde¬ 
pendent  on  the  inlet  temperature  T%,  until  a  maximum  value 
is  reached.  The  lower  the  lower  is  this  maximum  and  the 
lower  is  the  cell  voltage  E  where  the  maximum  is  obtained.  For 
cell  voltages  higher  than  the  maximum  value,  the  power  drops 


rapidly.  This  drop  is  again  almost  linear,  with  different  slopes  for 
different  temperatures.  This  result  supports  a  careful  choice  of 
the  operating  cell  voltage,  because  the  power  might  drop  by  more 
than  1  W  for  a  change  of  cell  voltage  of  only  0. 1  V,  especially 
for  higher  temperatures  T%.  In  Fig.  5(b),  the  exergetic  efficiency 
fi  is  given  as  a  function  of  the  cell  voltage  E  for  different  SOFC 
inlet  temperatures  7s- 

An  interesting  aspect  of  Fig.  5(a)  and  (b)  is  the  abrupt 
decrease  of  power  after  the  maximum,  giving  the  graphs  an 
almost  triangular  shape.  Fig.  5(c)  describes  one  reason  for  this 
phenomenon.  The  typical  voltage-current  curve  is  shown  for 
T%  =  800  K.  At  a  cell  voltage  E  of  about  0.55  V,  the  slope  of  the 
voltage-current  curve  changes.  This  happens  near  the  voltage 
where  the  maximum  power  P  is  achieved.  For  higher  voltages 
than  0.55  V,  the  curve  is  dominated  by  activation  overpotentials, 
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Fig.  6.  (a)  Electric  power  P  as  a  function  of  the  air/fuel  ratio  A.sofc  for  different 
SOFC  inlet  temperatures  7g.  Labels  indicate  T%  in  K.  All  other  operational 
parameters  are  taken  from  the  standard  case,  (b)  Exerg etic  efficiency  /i  as  a 
function  of  the  air/fuel  ratio  A.sofc  for  different  SOFC  inlet  temperatures  T%. 

whereas  for  smaller  cell  voltages,  the  curve  is  mainly  defined  by 
ohmic  losses.  This  effect  is  intensified  by  another  phenomenon, 
shown  in  Fig.  5(d).  This  plot  presents  the  current  density  profiles 
i  along  the  fuel  cell  for  Tg  =  800  K  and  different  cell  voltages  E. 
For  all  cell  voltages,  the  maximum  current  density  is  achieved 
close  to  the  beginning  of  the  fuel  cell.  For  smaller  voltages,  the 
current  density  drops  dramatically  to  zero  at  a  certain  point, 
because  all  hydrogen  is  consumed.  The  higher  the  cell  voltage, 
the  more  the  point  with  i  &  0  moves  towards  the  end  of  the  fuel 
cell.  For  cell  voltages  £>0.55  V,  the  fuel  cell  is  too  short  to 
convert  all  hydrogen. 

Another  important  operating  parameter  is  the  air/fuel  ratio 
Asofc.  In  Fig.  6(a),  the  electric  power  P  is  plotted  as  a  function 
of  the  air/fuel  ratio  A.sofc  for  different  SOFC  inlet  temperatures 
TV  The  labels  indicate  the  temperature  TV  The  electric  power 


Fig.  7.  (a)  Electric  power  P  as  a  function  of  the  heat-transfer  coefficients  a £  and 
ajj  for  different  SOFC  inlet  temperatures  7g.  Labels  indicate  7g  in  K.  All  other 
operational  parameters  are  taken  from  the  standard  case,  (b)  Exergetic  efficiency 
ji  as  a  function  of  the  heat- transfer  coefficients  a £  and  ajj  for  different  SOFC 
inlet  temperatures  T%. 

increases  linearly  if  Asofc  <0.5,  independent  on  temperature.  In 
this  region,  the  air  supply  determines  the  achievable  power.  If 
n  >  900  K,  a  maximum  power  of  about  5.5  W  is  reached  for 
air/fuel  ratios  slightly  higher  than  1.0.  Further  increase  of  the  air 
inlet  flow  rate  leads  to  no  further  increase  in  power,  since  the 
performance  is  limited  by  the  hydrogen  supply  in  this  operating 
region.  For  lower  temperatures  TV  the  air/fuel  ratio  affects  the 
performance  more  severely.  If  Ts  -  800  K,  a  maximum  power 
of  5.5  W  can  be  reached,  if  Asofc  is  increased  to  a  value  of  9.0. 
A  high  air  inlet  flow  rate  can  therefore  decrease  the  necessary 
SOFC  inlet  temperature  T%  significantly.  Fig.  6(b)  shows  the 
exergetic  efficiency  gasa  function  of  the  air/fuel  ratio  A.sofc  for 
different  SOFC  inlet  temperatures  TV 
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Fig.  8.  (a)  Electric  power  P  as  a  function  of  the  inlet  molar  flow  rate  of  butane 
^c4h10,i  f°r  different  SOFC  inlet  temperatures  T%.  Labels  indicate  T%  in  K.  All 
other  operational  parameters  are  taken  from  the  standard  case,  (b)  Exergetic 
efficiency  /jl  as  a  function  of  the  inlet  molar  flow  rate  of  butane  72c4Hi0,i  for 
different  SOFC  inlet  temperatures  7g. 

For  the  standard  case,  the  heat-transfer  coefficients  between 
environment  and  cathode  channel  a £  and  anode  channel  a £  are 
both  set  to  1  W  m~2  K_1 .  These  parameters  have  a  strong  influ¬ 
ence  on  the  electric  power  and  the  efficiency  of  the  fuel  cell 
system.  In  Fig.  7(a),  the  electric  power  P  is  plotted  as  a  function 
of  the  heat-transfer  coefficients  and  for  different  SOFC 
inlet  temperatures  T%.  The  labels  indicate  the  temperature  7% 
All  other  operational  parameters  are  taken  from  the  standard 
case  presented  in  Table  1 .  The  star  marks  the  result  of  the  stan¬ 
dard  case.  For  small  heat-transfer  coefficients,  indicating  a  good 
insulation  of  the  fuel  cell,  a  maximum  power  of  about  5.45  W 
is  reached  for  all  temperatures.  For  a  low  SOFC  inlet  tem¬ 
perature  of  800  K,  the  power  P  starts  to  drop  for  heat-transfer 
coefficients  below  0.5  W  m-2  K_1 .  If  T%  >  900  K,  the  drop  starts 


above  2.0  W  m-2  K_1 .  The  performance  loss  is  remarkable.  The 
achievable  power  P  can  almost  be  doubled  for  T%  =  800  K  by 
improving  the  insulation  from  1.0  to  0.5  W  m~2  K_1.  This  may 
offer  another  possibility  to  decrease  the  operating  temperature 
of  the  SOFC  without  losing  electric  power.  In  Fig.  7(b),  the  exer¬ 
getic  efficiency  /x  is  presented  as  a  function  of  the  heat-transfer 
coefficients  g£  and  a £  for  different  SOFC  inlet  temperatures  7s. 

A  final  analysis  shows  the  dependence  of  the  performance 
of  the  entire  fuel  cell  system  on  the  molar  flow  rate  of  liquid 
butane  ftc4Hi0,i-  In  Fig.  8(a),  the  electric  power  P  is  plotted  as 
a  function  of  the  inlet  molar  flow  rate  of  liquid  butane  ftc4Hi0,l 
for  different  SOFC  inlet  temperatures  7s.  The  labels  indicate 
the  temperature  7%  All  other  operational  parameters  are  taken 
from  the  standard  case  presented  in  Table  1.  The  star  marks 
the  result  of  the  standard  case.  For  SOFC  inlet  temperatures 
T 8  >  900  K,  the  power  increases  almost  linearly  with  the  molar 
flow  rate  of  butane,  since  all  hydrogen  is  converted  in  the  fuel 
cell.  For  lower  temperatures,  the  performance  of  the  fuel  cell  is 
too  low  to  use  all  available  hydrogen.  In  Fig.  8(b),  the  exergetic 
efficiency  /x  is  shown  as  a  function  of  the  inlet  molar  flow  rate 
of  butane  /?c4Hi0,i  for  different  SOFC  inlet  temperatures  T%.  For 
T%  =  900  K,  the  exergetic  efficiency  /x  is  almost  constant  within 
the  plotted  range. 

The  results  of  this  study  show  that  the  fuel  cell  performance 
can  be  easily  increased  by  adjusting  the  operating  parameters. 
For  the  standard  case  with  a  SOFC  inlet  temperature  T%  =  900  K, 
the  power  P  can  be  enhanced  by  slightly  changing  the  cell  volt¬ 
age  E  to  0.78  V  and  the  air/fuel  ratio  Asofc  to  1.2.  These  small 
adjustments  lead  to  an  optimum  electric  power  of  6. 1  W,  com¬ 
pared  to  5.45  W  in  the  standard  case.  The  exergetic  efficiency  /x 
rises  to  21.9%.  The  results  demonstrate  possibilities  to  achieve 
satisfactory  power  for  SOFC  inlet  temperatures  T%  <  900  K.  By 
using  the  parameters  of  the  standard  case  and  only  changing  the 
SOFC  inlet  temperature  T%  to  800  K,  the  power  drops  to  2.9  W. 
However,  by  changing  the  cell  voltage  to  0.55  V,  increasing  the 
air/fuel  ratio  Asofc  to  2.0,  and  improving  the  heat-transfer  coeffi¬ 
cients  and  ajfJ  to  0.5  W  m-2  K-1,  the  power  P  reaches  4.3  W. 
The  exergetic  efficiency  amounts  to  15.4%  in  this  case. 

5.  Conclusions 

The  fuel  cell  micropowerplant  model  developed  in  this  study 
replicates  the  characteristic  behavior  of  a  SOFC  system  under  a 
host  of  operating  conditions.  The  influence  of  operational  and 
geometric  parameters  on  the  efficiency  of  a  micro  SOFC  sys¬ 
tem  is  investigated  according  to  the  first  and  the  second  law 
of  thermodynamics.  It  is  shown  that  the  exergetic  efficiency 
of  the  system  can  be  increased  significantly  by  proper  choice 
of  operating  conditions.  Useful  information  about  such  opti¬ 
mal  conditions  can  be  deduced  from  the  model  presented  in 
this  study.  The  developed  fuel  cell  model  can  be  modified  and 
extended  to  consider  different  types  of  fuels  or  different  sys¬ 
tem  configurations.  Further  improvement  could  be  achieved  by 
extending  the  simplified  flow  model  in  the  fuel  cell  channels 
and  by  using  a  more  complex  heat- transfer  model.  Neverthe¬ 
less  all  relevant  physics  is  covered  by  the  present  model.  An 
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interesting  approach  might  be  the  simultaneous  conversion  of 
different  species  in  the  fuel  cell  instead  of  single  usage  of  hydro¬ 
gen. 

Appendix  A 


hydrogen  and  water,  given  as  [9] : 


f^eff ,  a  — 


cH20,h^eff,H2  +  cH2,h^eff,H20 


.a 


.a 


cH20,h  +  cH2,h 


The  diffusion  coefficient  of  hydrogen  [9], 


(A. 7) 


The  activation  polarizations  ^act,a  and  r]act,c  are  expressed  by 
Butler- Volmer  equations  [7,9,18].  The  activation  polarization 
of  the  anode  and  cathode  is  given  as  [9] : 


7? act,  a 


2  RT 
n&F 


ar  sinh 


and 


7?act,c 


2  RT 
nQF 


ar  sinh 


(A.l) 


(A.  2) 


respectively,  where  ne  indicates  the  number  of  electrons  trans¬ 
ferred  in  a  single  elementary  rate-limiting  step  and  a  common 
assumption  is  zze  =  1  [8].  The  current  density  is  represented  by  i 
and  the  exchange  current  densities  z'o,a  and  z'o,c  are  the  electrode 
reaction  rates  at  the  equilibrium  potential.  In  order  to  account  for 
the  temperature  dependence  of  the  exchange  current  densities, 
the  exchange  current  densities  are  modeled  as  [14]: 


,-o 

0,a 

T 


exp 


for  the  anode,  and 


,-o 

0,c 

T 


(A. 3) 


(A.4) 


£>eff,H2  —  (^eff,H2,^  +  ^eff,  H2-H2C>)  (A.8) 

and  the  diffusion  coefficient  of  water  [9] 

£>eff,H20  =  (^eff,H20,K  +  ^eff,H2-H2o)  (^.9) 

represent  a  combination  of  Knudsen  diffusion  and  binary  ordi¬ 
nary  diffusion.  The  Knudsen  diffusion  coefficients  [19], 


2  i cy  8 RT 

%,h2,k=-x10t-^—  (A' 10) 

and 


2  ,  c  y  %RT 

=  j  X  10  ■  (A.  II) 

of  hydrogen  and  water,  respectively,  are  functions  of  the  tem¬ 
perature  T,  the  molar  mass  of  the  gas  M,  and  the  average  pore 
radius  r  of  the  anode.  In  order  to  take  into  account  the  tortuous 
path  of  the  diffusing  gases  and  the  porosity  of  the  electrodes,  the 
effective  diffusion  coefficients  are  multiplied  with  the  ratio  of 
the  porosity  y  and  the  tortuosity  §.  The  binary  ordinary  diffusion 
coefficient  in  m2  s-1  is  calculated  according  to  the  Chapman- 
Enskog  equation,  given  as  [20] : 


for  the  cathode,  using  Arrhenius  type  coefficients  /q  a,  /[ j  c,  £a, 
and  Ec  extracted  from  Table  1  of  Ref.  [14]. 

The  concentration  polarizations  77Conc,a  and  r]conc,c  describe 
the  resistances  for  mass  diffusion  through  the  porous  electrodes. 
In  the  anode,  ID  equimolar  diffusion  of  hydrogen  and  water 
takes  place  at  steady  state,  where  the  molar  fluxes  are  functions 
of  molar  concentration  gradients  and  the  effective  binary  dif¬ 
fusion  coefficient  DQ ff?a.  By  combining  the  definitions  of  molar 
fluxes  of  hydrogen  and  water,  the  concentration  polarization  for 
the  anode  is  [9] : 


7?conc,a  — 


RT 

~1F 


In 


(  1  -  (4/2FDeff,ac4  h)*  \ 

\  1  +  (4/2 FDeff.aC^o^)/ )  ’ 


(A. 5) 


A:ff,H2-H2o  =  1-8583  x  10  7 


rl5([l/MH2]  +  [l/MH2o])1/2 
^aH2— h2o  ^2h2-H2  o 


(A.  12) 


Herein,  p  is  the  pressure  in  atmospheres.  The  collision  diam¬ 
eter  cth2-h2o  with  units  in  angstrom  and  the  collision  integral 
o  are  based  on  the  Lennard- Jones  potential  [21].  Simi¬ 
larly,  the  effective  diffusion  coefficient  of  the  cathode  Deff,c  is  a 
combination  of  Knudsen  diffusion  and  binary  ordinary  diffusion 

[9] 


Ajff,c  =  7Jeff,02  =  (D 


-l 

eff,02,K 


+  D 


-l 

eff,02— N2 


(A.  13) 


where  /a  is  the  anode  thickness.  The  concentration  polarization  and  8o2  is  the  ratio  of  these  two  diffusion  parts  [9]: 
of  the  cathode  reads  [9] : 


RT 

7?conc,c  — 


(p/8q2)  -  (( p/8q2 )  -  RTcc0i  h)  exp(RT8o2lci/4FDQ^cp) 


RT<-0^ 


(A.6) 


where  lc  is  the  cathode  thickness  and  Cq2  h  the  molar  concentra¬ 
tion  of  oxygen  in  the  cathode  channel.  The  factor  5o2  is  the  ratio 
of  Knudsen  diffusion  and  the  sum  of  binary  ordinary  diffusion 
and  Knudsen  diffusion.  The  effective  diffusion  coefficient  of  the 
anode  £>eff?a  is  calculated  by  using  the  diffusion  coefficients  of 


A:ff,02,K 

£*eff,02,K  +  £)eff,o2-N2 


(A.  14) 
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The  Knudsen  diffusion  coefficient  of  oxygen  is  given  as  [19]: 


£>eff,02,K 
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7TMq2 


(A.  15) 


and  the  binary  ordinary  diffusion  coefficient  in  m2  s  1  reads 

[20]: 


Ajff,o2-N2  =  1.8583  x  10  7 


rl5([l/Mo2]  +  [l/MN2])1/2 

Pa O2-N2 


X 


(A.  16) 


The  last  polarization  considers  ohmic  losses  due  to  transport  of 
ions  in  the  electrolyte.  Ohmic  losses  in  the  electrodes  are  much 
smaller,  and  are  therefore  neglected.  The  electric  current  in  the 
electrolyte  obeys  Ohm’s  law, 


7? ohm  —  RqP- 


(A.  17) 


Taking  into  account  the  dependence  on  temperature,  the  specific 
ohmic  resistance  of  the  electrolyte  RQ\  can  be  expressed  as  [7]: 


^el 

<70 


(A.  18) 


where  /e  1  refers  to  the  thickness  of  the  electrolyte,  Ee\  the  acti¬ 
vation  energy  of  the  ionic  transport,  and  a 0  a  reference  ionic 
conductivity.  The  coefficients  ao  and  EQ\  are  extracted  from  Fig. 
8  of  Ref.  [15]. 
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